FERMILAB-Pub-93/399-A 
SUSSEX-AST 94/2-1 
astro-ph/9402021 
(February 1994) 
submitted to Physical Review D 



SECOND-ORDER RECONSTRUCTION 
OF THE INFLATIONARY POTENTIAL 

Andrew R. Liddle 1 and Michael S. Turner 2 ' 3 

1 Astronomy Centre, Division of Physics and Astronomy 
University of Sussex, Brighton BN1 9QH, UK 

2 NASA/Fermilab Astrophysics Center 
Fermi National Accelerator Laboratory, Batavia, IL 60510-0500 

3 Departments of Physics and of Astronomy & Astrophysics 
Enrico Fermi Institute, The University of Chicago, Chicago, IL 60637-1433 

ABSTRACT 

To first order in the deviation from scale invariance the inflationary potential and its first 
two derivatives can be expressed in terms of the spectral indices of the scalar and tensor 
perturbations, n and n-r, and their contributions to the variance of the quadrupole CBR 
temperature anisotropy, S and T. In addition, there is a "consistency relation" between these 
quantities: rix = — 1§ • We discuss the overall strategy of perturbative reconstruction and 
derive the second-order expressions for the inflationary potential and its first two derivatives 
and the first-order expression for its third derivative, all in terms of n, tit, S, T, and dn/dln k. 
We also obtain the second-order consistency relation, n T = — =^[1 + 0.11^ + 0.15(n — 1)]. 
As an example we consider the exponential potential, the only known case where exact 
analytic solutions for the perturbation spectra exist. We reconstruct the potential via Taylor 
expansion (with coefficients calculated at both first and second order), and introduce the 
Pade approximant as a greatly improved alternative. 



PACS number(s): 98.80. Cq, 98.70.Vc 



1 Introduction 



In inflationary models quantum fluctuations excited on very small length scales (~ H _1 ~ 
10~ 23 cm) are stretched to astrophysical scales (~ 10 25 cm) by the tremendous growth of the 
scale factor during inflation (H is the value of Hubble parameter during inflation) 0] . This 
results in almost scale-invariant spectra of scalar (density) and tensor (gravitational- 
wave) metric perturbations. Together with the prediction of a spatially-flat Universe 
they provide the means for testing the inflationary paradigm. The tensor fluctuations lead to 
cosmic background radiation (CBR) anisotropy and a stochastic background of gravitational 
waves with wavelengths from about 1 km to over 10 4 Mpc. The scalar fluctuations also lead 
to CBR anisotropy and seed the formation of structure in the Universe. 

The amplitudes and spectral indices of the metric fluctuations can be expressed in terms 
of the inflationary potential and its derivatives, evaluated at the value of the scalar field when 
astrophysically interesting scales crossed outside the horizon during inflation (from galactic 
scales to the presently observable Universe, corresponding to the eight e-foldings about 50 
e-folds or so before the end of inflation). Techniques have been developed for relating the 
scalar and tensor spectra to the potential and its derivatives in an expansion whose small 
parameter is the deviation from scale invariance |4], ||. In particular, the spectral indices 
and the power spectra of the fluctuations today can be written as 
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Here k is the comoving wavenumber, x = mp\V'/V measures the steepness of the potential, 
prime denotes derivative with respect to the scalar field that drives inflation, subscript 50 
indicates that the quantity is to be evaluated 50 e-folds before the end of inflationj 2 ] mp\ = 

1 Several minor errors in Ref. § have been corrected here: the factor of H^ +n in Eq. (A5) should be Hq ; 
the factor of H^ +n in Eq. (A7) should be 2™ _1 7Jq; the factor of l.l(n — 1) in Eq. (A8) is more precisely 
1.3(n — 1); the factor of 1.2tit in Eq. (A14) is more precisely \Aut- 

2 The point about which the potential is expanded is in principle arbitrary. However, the spectral indices 
n and tlt can only plausibly be measured on scales from 1 Mpc — 10 4 Mpc and S and T depend upon 
perturbations on these same scales, so it makes sense to choose the expansion point to correspond to when 
these scales crossed outside the horizon during inflation; in addition, by taking k 50 T = 1 several expressions 
simplify. The precise number of e-folds before the end of inflation when these scales crossed outside the 
horizon depends logarithmically upon the energy scale of inflation and the reheat temperature, see Refs. §, 
H ^J; for the sake of definiteness we take this number to be 50, which can easily be changed to the correct 
value for a given inflationary model. 
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1.22 x 10 19 GeV is the Planck mass, Hq is the present value of the Hubble constant, tq ~ 2if ~ 1 
is the present conformal age of the Universe, and 7 ~ 0.577 is Euler's constant. Scale- 
invariant metric perturbations correspond to (n — 1) = ut = 0. The functions T(k) and 
Tx(k) are the transfer functions for scalar [|7| and tensor || metric perturbations respectively; 
for kro <C 100, both T(k) and Tx{k) — » 1. The expressions for n and tit are given to 
lowest order in the deviation from scale invariance (hereafter, referred to as first order), and 
the expressions for A and At include the lowest-order term as well as the next correction 
(hereafter, referred to second order) [[|. 

From these expressions the consequences of the scalar and tensor metric fluctuations may 
be computed. In particular, the contributions to the variance of the angular power spectrum 
of the CBR anisotropy on large angular scales (I <C 200) which arise predominantly due to 
the Sachs-Wolfe effect are given by || 

(l«fj 2 > = ir k- 2 p(k)\MkT )\ 2 dk, (5) 
(l«Ll 2 ) = ae^^^ATjT^i^wHTrC*)! 2 ^, (6) 

where t lss ~ t /(1 + z LSS ) 1//2 ~ t /35 is the conformal age at last scattering (zlss — 1100) 
and ji is the spherical Bessel function of order I. (We note in passing that both expressions 
are based upon the approximation that the Universe is matter-dominated at last-scattering; 
the small contribution of radiation, about 10%-20%, leads to corrections || that would have 
to be included in a more accurate treatment. The corrections to the quadrupole anisotropy 
are small.) 

The contribution of scalar and tensor metric perturbations to the observer averaged 
variance of the quadrupole CBR anisotropy can be computed numerically || 

S = 5 -Ml£l ~ 2.2 [1 + 1.2n T + 0.08(n - 1)] V ^ k ff~\ (8) 

47T ^P1 X 50 

t . 5 MI1 ^ 0Mll + 1AnT] M^ri t (9) 
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where the dependence upon (n — 1) and rir is given to first-order. In evaluating these 
expressions the effect of transfer functions is negligible as the integrals are dominated by 
kro ~ 2. For simplicity, following footnote 2 we henceforth omit factors of (fcsoTo) 1- ™ and 
(k 50 T )~ nT ; they are easily re- inserted if needed. 
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1.1 First-order reconstruction 



We choose S, T, n T , and (n — 1) as a convenient set of observables; other choices are possible 
and can be easily transformed to our set.Q Since S, T, n, and tit are expressed in terms of the 
potential and its first two derivatives, one can invert the expressions to solve for the potential 
and its first two derivatives in terms of S, T, n, and ut plus a "consistency relation." Those 
expressions are [[12] 
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In the second expressions for the potential and its first two derivatives we have used the 
consistency relation to express ut in terms of as ^ should be easier to measure [|10j. 
Note that the sign of V cannot be determined as it can be changed by a field redefinition 
— > —0, though a specific choice here determines the signs of various later expressions. This 
procedure actually generates the full second-order term for V50, while the other expressions 
are first-order. 

In order to actually reconstruct the inflationary potential over the eight or so e-folds 
relevant for astrophysics from its value and first two derivatives one needs to relate the 
number of e-foldings from the end of inflation, N (where dN/dt = —H), to the value of the 
scalar field. To lowest order the equation for d<fi/dN follows from the slow-roll equation for 
the evolution of <t> and is given by 
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where to lowest order the right-hand side is just yj—nr mp\/ V8vr. 

In the next section we discuss the overall strategy of perturbative reconstruction, and in 
the following section go on to derive the second-order expressions for the potential and its 
first two derivatives, for the equation relating and N, and for the consistency relation, as 
well as the first-order expression for the third derivative. We finish with a brief discussion 
of our results. 



3 For example, in separating the tensor and scalar contributions to CBR anisotropy one might measure 
/(/+l)(|a/ m | 2 ) for four values of I (or ranges centered on four different values of I); see e.g., Refs. JlO|, fl). From 
these measurements and the known dependence of 1(1 + l)(|a^ 



T{l/2) nT upon S, (n - 1), T, and n T , cf. Eqs. (g) and 



)/4tt«5(Z/2)"- 1 and 1(1 + l)(\afj 2 ) /4n 
our chosen observables can be extracted. 
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2 Perturbative Reconstruction Strategy 



While one can hope to learn about the potential over the interval that affects astrophysical 
scales, it is probably not realistic to hope to learn much about the potential globally without 
some additional a priori knowledge (e.g., the functional form of the potential). The funda- 
mental goal of perturbative reconstruction is to use a finite set of data to reconstruct the 
inflationary potential over the interval where the eight or so e-foldings of inflation relevant 
to astrophysics took place. The observational data all trace to the scalar and tensor metric 
perturbations, whose observable consequences can be expressed in terms of the inflationary 
potential V and its derivatives evaluated at some convenient point in this interval. (For 
brevity, in this section we drop the subscript '50' that indicates where the potential and its 
derivatives are to be evaluated.) Once the observables, e.g., n T , (n — 1), S, T and so on, 
have been expressed in terms of the potential and its derivatives, these expressions can be 
inverted to express the potential and its derivatives in terms of the observables, as well as a 
consistency relation. From these the potential can be recovered by expansion. 

In principle, the observables depend upon all the derivatives of the potential, making the 
problem appear intractable. If one is willing to restrict the problem to flat potentials which 
lead to nearly scale-invariant perturbations and nearly exponential inflation, the problem can 
be made manageable. (In the scale-invariant limit the potential is precisely constant and 
all its derivatives vanish.) In the nearly scale-invariant limit we have a set, albeit infinite, 
of small parameters to expand in: m^y'™' /V; as we shall describe, when calculating to a 
given accuracy only a small number of derivatives are needed. Put another way, the terms 
involving more derivatives or higher derivatives are of higher order. To be more specific, to 
lowest order 

T ~ 0(V/m Pl ), (15) 
n T , T/S ~ 0[(m P1 V'/V)% (16) 
(n-1) ~ 0[(m Pl V'/V) 2 ] + 0[m 2 Pl V"/V}; (17) 

that is, only the potential and its first two derivatives come into play. 

Corrections from higher derivatives come into play because of the variation of the poten- 
tial and its first two derivatives during the Hubble time or so that a given scale is crossing 
outside the horizon and is becoming a classical metric perturbation. It is straightforward 
to write down the form of the higher derivative terms expected by using the fact that the 
variation of a given derivative over a Hubble due to a higher derivative is: 

SV (n) ~ V {m) 6(f) m - n ~ {ml x V' /V) m - n V {m \ (18) 

where the final expression follows by using d<p/dN ~ mpiV'/V. The form of the higher-order 

4 The one possible exception involves the accurate measurement of the stochastic background of gravita- 
tional waves on scales from 1 km to 3000 Mpc (corresponding to N ~ — 50) in which case the inflationary 
potential could be mapped out directly since the amplitude of the tensor perturbation on a given scale is 
related to the value of the potential. 
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terms in the expansions of T, ut, T/S, and n can now be written down directly: 

T ~ 0(V/m Pl )[l + 0[(m Pl V'/V) 2 ] + 0[(m 2 Pl V''/V)(m Pl V'/V) 2 ] + --], (19) 

n T ,T/S ~ ©[(mp^'/yj^+o^yvvoCmpiVV^) 2 ] 

+ 0[(m P1 \/( 3 VV r )(m P1 y , /y) 3 ] + ---, (20) 
n - 1 ~ ©[(mpiF'/y) 2 ] + 0[m P1 F"/y] + 0[(m P1 V (3) /VO(mpiV7V)] 

+ 0[(m Pl V®/V)(m Pl V"/V)(m P1 V'/V)] + 0[(m P1 ^ 3 )/y)(m P1 y7F) 3 ] + • (21) 

The expansion for V begins with a term that involves no derivatives; the next term 
involves two derivatives; the next four derivatives, and so on. The expansions for nj>, T/S, 
and (n — 1) begin with terms involving two derivatives; followed by terms involving four 
derivatives; and so on. In the previous literature, the lowest-order term has been referred 
to as first-order; the next term, which involves two additional derivatives, has been referred 
to as second-order; and so on. Explicit expressions for the second-order terms are given in 
Refs. 0, T^]; some of the third-order terms for (n — 1) are given in Ref. 



In the next three subsections we address the convergence of the Taylor series for the 
potential and the relative sizes of the terms in the expansions for the observables. We 
show that for a very general class of potentials that lend themselves to reconstruction that 
the higher-order terms in these expansions are smaller and are bounded by nx/AiV™/ 2-1 , 
where m is the number of derivatives in the term, and further, that the Taylor series for the 
potential is absolutely convergent. 

Before going on, let us remind the reader of a very useful fact and mention some notation. 
The variation in the scalar field over the eight relevant e-folds of inflation will be needed in 
many places; it is A0/m P i ~ (m-piV /V)AN, where AiV ~ 8 and throughout we use A to 
indicate the change in a quantity over the eight relevant e-folds. Since tit ~ (m P iV /V) 2 , 

1 /2 

we will use n T to characterize the size of m P {V /V . While it is actually (n — 1) — 3n T and 
not {n — 1) whose lowest order term is given by m Pl V"/V, for simplicity we will often use 
(n — 1) to characterize the size of m Pl V"/V. 

2.1 Scale-free potentials 

Let us begin with a very simple class of potentials before we consider the general case. These 
are potentials without a scale other than an overall normalization; e.g., V(4>) = Vq exp(— ft<f>), 
V(4>) = a<p b , or V{4>) = a<f)~ b . For such potentials there is but a single expansion parameter 
since 

my(m) 

0[(m Pl V'/V) m }. (22) 

For the potentials given above, Eq. (p2|) follows directly; in the absence of a more quantitative 
definition of scale-free we shall use Eq. ([23) as the definition. 

If we use ~ (m P \V' /V) 2 to characterize the deviation from scale-invariance, it follows 
that 

m P\ V «r m/2 



V 



0\n™'\ (23) 
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with higher-order corrections to the expression for mfy^ /V increasing as powers of tit- 
For such potentials n — 1 and are necessarily of the same order, which is not true in the 
general case. The convergence of the power series for V(4>) over the interval of AN ~ 8 
e-folds is manifest, as the contribution of the m-th derivative to the Taylor expansion is 

AV m 1 V^Acf)" 1 {n T AN) m 



V ml V ml 



For scale-free potentials Eq. (p^) provides the ordering of terms in the expansion of the 
observables in terms of the derivatives of the potential very directly: the order of a term 
involving m derivatives is (mp\V /V) m ~ 0(n™^ 2 )- For example, in expression Eq. for 
(n — 1) the first two terms are of the order of n^; the next is of the order of n^; and the final 



two are of the order of n\. 



2.2 Strong reconstructability 

A priori we do not know the form of the potential and thus whether or not it is scale- free; 
therefore, it is important to address the most generic case possible. Lacking a priori knowl- 
edge of the potential, one can take advantage of the observational data itself for guidance in 
reconstruction. In the near term the observational data available are likely to be a handful 
of numbers, e.g., tit, (n — 1), S, and T. A reasonable, robust, and pragmatic criterion for 
reconstructability is that the spectral indices do not vary greatly over the eight e-folds of 
interest; that is, 

|Anr| |A(n-l)| 
\n T \ \{n — l)\ 

where 5 is some suitably small number. We shall refer to this as "strong reconstructability," 
or SR. 

Since the scalar and tensor spectral indices depend upon the first two derivatives of the 
inflationary potential, SR can be quantified in the following way: V and V" should not vary 
significantly over the interval of inflation affecting astrophysically interesting scales. This 
in turn constrains the higher derivatives of the potential through their contributions to the 
Taylor expansions of V and V": 

AV V^Ad)" 1 - 1 V 

AV" V^ m) Ad) m - 2 V" , 

Again using A<p/m-p\ ~ (mpiV /V)AN ', these bounds become 



m 



V \ V 

< 0{n T m/2+1 AN- m+1 5{m-\y] for m > 2; (2* 
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m ^ V "\ AN~ m+2 5 (m-2)! 



m+2 / 2 

< 0[(n- l)nr m/2+1 AA^- m+2 5(m-2)!] for m > 3, (29) 



where in the final expressions we have used the fact that ny ~ (^[(mpiK'/V) 2 ] and that 
(n-1) ~0[m P1 V"/V). 

These constraints to the derivatives of the potential are weaker than the ones we derived 
for scale-free potentials, but are more generally applicable and serve the same purpose. The 
second of these implies that the Taylor expansion for the potential is absolutely convergent, 
as it bounds the contribution to V{<p) from the m-th derivative 

AV m V^A<f) m f n T (n-l)AN 2 5 \ 

~~T? 77} < U 7 7\ • (30) 

V mlv V m[m — 1) I 



Constraints (28) and (p9|) also serve to order terms in the expansions of the observables 
in terms of the potential and its derivatives. For example, Eq. ( ^IP for the scalar index 
includes terms of order {m P yV /V) 2 , m P iV"/V, {m P yV' /V)(m Pl V^ /V). Our SR bounds tell 
us nothing about the relative sizes of first and second terms (the first-order terms), though 
they imply that the second term must be smaller than 5 /AN. Based on the SR bounds, 
the third term (second-order term) must be less than both 5/ AN 2 and (m Pl V" /V)S / AN , 
and so it is necessarily of higher order than the second term. The higher- derivative terms, 
(m Pl V'/V) 3 (m Pl V^/V) and {m Pl V /V)(m Pl V" '/V)(rr4iV^ /V) terms are of even higher 
order: the former must be less than 5/ AN 2 times the first term in the expansion for (n — 1) 
and less than htS/AN times the second term in the expansion for (n — 1), while the latter 
must be less than 6/ 'AN 2 times the second term in the expansion for (n — 1). The ordering 
of the terms in the derivative expansion for (n — 1) is clear: more derivatives are suppressed 
by powers of AiV. In particular, a term involving m derivatives can be no larger than 
nrd/AN™/ 2 - 1 . 

Before ending this subsection, we mention an interesting possibility: for models with very 
large deviations from scale-invariance the data may some day be good enough that a small 
fractional change in the spectral index is observable. A case in point is intermediate inflation 
|15| , where the scalar index may be greater than unity and may decrease significantly. In 



particular, the potential for intermediate inflation is scale-free and dn/dlnk ~ (n — l) 2 so 
that A(n — l)/(n — 1) ~ AN(n — 1). The SR bounds still apply, and additionally, the new 
observable dn/d\nk allows one to determine V" at lowest order (as described in section 3). 



2.3 Weak reconstructability 

The pragmatic criteria of SR discussed above can be relaxed somewhat, without sacrificing 
the convergence of the Taylor series for the potential or the ordering of terms in the expan- 
sions for the observables. Suppose that one, or even both, of the spectral indices did indeed 
exhibit a large fractional change over astrophysically interesting scales, so that a power- 
law description of the scale dependence of the metric perturbations is not strictly valid. If 
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the absolute value of the change is much less than unity, then the fact that Aut > ut 
or A(n - 1) > (n — 1) is undetectable and of little practical significance, and, as we shall 
show now, reconstruction can proceed. We refer to this as "weak reconstructability," or 
WR. A case in point is the natural inflation model [1£|; with parameters chosen to give 



(n — 1) = —0.3, the tensor spectral index grows by a factor of about 100 between the largest 
and smallest interesting scales. However, this growth is entirely unobservable, being the 
difference between rij- = — 10 -9 and = — 10~ 7 . 

Logically, there are three cases of WR: (i) scalar index satisfies SR and tensor index 
satisfies WR; (ii) tensor index satisfies SR and scalar index satisfies WR; and (iii) both tensor 
and scalar indices satisfy WR. Since we have previously derived the bounds to rripyV^ /V 
that follow if tensor and scalar indices satisfy SR, cf. Eqs. (|28D and (|29|) respectively, here 
we simply do the same for WR. In case (i) the SR scalar and WR tensor bounds apply; in 
case (ii) the WR scalar and SR tensor bounds apply; and in case (iii) the scalar and tensor 
WR bounds apply. In all three cases the implications for convergence of the power series 
and the ordering of terms is very similar to the case of SR. 

Let us take S to be the parameter that quantifies the smallness of the tensor (or scalar) 
index and its absolute change. For sake of definiteness, we would imagine that a change of 
a few hundredths for the scalars, and considerably more for the tensors, would be extremely 
hard to observe. Following the same strategy as in the SR case, this time bounding the 
absolute change in the spectral index due to higher-order derivatives, we find: 



~v 



Pl V \ ^ „~ m / 2 A AT-m+1 



< rij "AN- m+l 5(m - 1)! for m > 2 (tensor), (31) 



j < n ( T- m)/2 AN- m+2 5{m - 2)! for m > 3 (scalar). (32) 

These constraints differ from their counterparts in the SR case only slightly: by one fewer 
factor of ft? (tensor) and by the absence of the (n — 1) factor (scalar). Thus, the conclu- 
sions reached for convergence and term ordering in the SR case carry over with only minor 
modification. For example, the size of the contribution of the m-th derivative to the Taylor 
series of the potential is bounded by (n — l)n^AiV 2 5/m(m — 1), nxAN 2 5 jmim — 1), and 
nTAN 2 5/m{m — 1) in cases (i)-(iii) respectively; this guarantees absolute convergence. If 
the SR tensor bound applies then a term involving m derivatives is as before bounded by 
utS j ' AN m l 2 ~ l ] if the WR tensor bound applies then such a term is bounded by 5/ AN" 1 ^ 2 ^ 1 . 

Finally, what types of potentials give rise to order unity fractional changes in the spectral 
indices while still satisfying the WR criteria? It is simple to show for the tensor index that 
AN(n — 1) must be of order unity or larger; this occurs in the previously mentioned natural 
inflation model. For the scalar index the condition is that (m^V^ /V)utAN / (n — 1) must 
be of order unity or larger. 
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2.4 Sensibility summary 

When physicists construct an expansion in a small parameter (or even several small param- 
eters) they rarely worry about rigorous mathematical issues. While we would like to follow 
in that tradition, the problem here is a bit more vexing as there are in principle an infinity 
of small expansion parameters: m^V^ /V. We have addressed two (not unrelated) issues 
here: convergence of the reconstructed potential and ordering of terms. 

Based upon pragmatic criteria that derive from the data themselves we have shown that 
convergence and term ordering follow for potentials where the spectral indices do not vary 
significantly over astrophysically interesting scales (referred to as SR), or if they do vary by 
of order unity, the absolute change is small by comparison to what can be measured (referred 
to as WR). In both cases we explicitly showed that the Taylor expansion for the potential is 
necessarily convergent, and that higher- derivative terms in the expansions for the observables 
descend in size. For "scale-free" potentials a term that involves m derivatives is of the order 
of n™^ 2 ; in the more generic cases of SR and WR, such a term is bounded by titS/AN 171 ^ 2 ^ 1 
and 5 / AN m / 2 ~ l respectively. This establishes what has been previously assumed implicitly 
in Refs. ||, |13| : the terms involving more derivatives are of higher order. 

2.5 The consistency relation 

An important feature of reconstruction is that the problem is overdetermined; specifically, 
a set of M > 3 observables can be expressed in terms of the potential and its first M — 2 
derivatives. This implies a "consistency relation," which, for increasing M, contains terms 
of higher and higher order. The lowest-order consistency equation, tit = — has been 
much discussed (e.g., in Ref. ||, ||]) and arises through Eqs. (|l|), (§) and (|) which express 
n T , S, and T in terms of V50 and V^. 

Calculating higher derivatives alone, while keeping the calculation of each derivative to 
lowest order, does not lead to the correct second-order term in the consistency equation, 
and nor does calculating the second-order corrections to the derivatives present. One must 
systematically do both. The second-order version of the consistency equation is obtained by 
calculating the potential, its derivative and Eq. (|I|) to a higher order. Adding an extra order 
to the calculation of V^ adds a new observable, {n — 1), which will appear in the consistency 
equation at second order. To account for there being still only a single consistency equation, 
there must be a new equation, and because (n — 1) has only entered at second-order in V^q, 
we only need the first-order equation for V^ . The second-order consistency equation, which 
we calculate in this paper, therefore relates ^ and (n — 1), with the last only appearing 
as a second-order correction. Were one to desire a calculation to yet higher order, the same 
pattern would persist; each existing derivative must be calculated to one extra order and 
the next derivative to lowest order, introducing a new observable. This will generate next- 
order terms in the consistency equation with the new observable appearing at that order. 
However, this presently cannot be done as third-order expressions for V50 and V 5 ' have not 
been calculated. 
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2.6 Expansion techniques 



Given the value of the potential and its first two or three derivatives at a point and the 
4>n relation just obtained, one can reconstruct the potential on the observationally relevant 
scales (i.e., N ~ 42 — 50). The standard technique used previously is the Taylor expansion 

V ((j)) = V 50 + V 5 ' o (0 - 5O ) + \VM - 5O ) 2 + • ■ ■ (33) 



For many situations this is perfectly fine (e.g., when ny and n — 1 are small, see Ref. ||12|| ). 
However, if the range of eight or so e-foldings corresponds to a large range in the conver- 
gence may not be very good because of the abrupt truncation of the Taylor series. Specifically, 
for large (<f) — 0so) the shape of the reconstructed potential is dictated, rightly or wrongly, 
by the last term in the expansion (quadratic or cubic). 

An alternative is the Pade approximant fTTH , which can be generated directly from a 
truncated power series. For a power series that extends to order N, the Pade approximants 
are quotients of two polynomials of order L (numerator) and M (denominator) denoted by 
[L, M], where L + M = N. By construction, the expansion of [L, M] matches that of the 
power series to order N, but of course is not truncated. Very often, Pade approximants 
provide a very good approximation over a wider range of values than the Taylor series from 
which they are derived; they in some way encode better estimates of the higher-order terms 
than does truncation. If we truncate the Taylor series at the second derivative, then the 
associated diagonal Pade approximant [1, 1] is a ratio of two first-order polynomials given 

byQ 

„,,v a o + ai(0-05o) , , 

1 + &l(0 - 050 ) 

with 

a = V 50 ; h = -V^/2V^ a, = - V 50 V>' /2V> . (35) 

As we shall illustrate later by specific example, Pade approximants have a lot to offer when 
the Taylor series proves a poor approximation. 



3 Second-order Reconstruction Reduced to Practice 

Having discussed the philosophy and strategy, let us proceed to deriving the full reconstruc- 
tion equations at second-order. The reconstruction equations for the scalar potential and its 
first two derivatives, evaluated to second-order, are given in Ref. [|13|], though not in terms 
of cosmological observables. They are given in terms of the perturbation amplitudes Aq 
and A 2 S . Very roughly, As is the horizon-crossing amplitude of the density perturbation on 
a given scale and Ac is the horizon-crossing amplitude of the tensor perturbation (in the 

J The [2, 0] approximant is just the truncated Taylor series; in addition to simplicity, there is some motiva- 
tion for using the diagonal approximant rather than the [0, 2] approximant as it is asymptotically constant, 
consistent with the flatness of inflationary potentials. 
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Appendix we provide some relations between notation used in that paper and this one.) Our 
purpose here is to express these second-order expressions for the potential and its first two 
derivatives in terms of the measurable quantities n, dn/dlnk, tit, S, and T. 
The amplitudes A 2 S and Aq are related to the observables S, T, tit and n by: 

Aq = 0.70(1 — 1.3ny)T, A 2 S = 9.6[1 - 1.15(n - 1)}S, (36) 

where the (n — 1) and dependencies have been found by evaluating the Sachs- Wolfe 
integrals numerically. Both expressions are accurate to second-order. 

Before deriving second-order expressions for the potential and its derivatives, we calculate 
the second-order version of the consistency relation. It is obtained from Eq. (2.9) of Ref. ||13|| , 

~Y = § S [ 1 + 3e - 2 ^ ( 3? ) 

where to the required order the slow-roll parameters e and rj (defined in the Appendix) are 
given by 

e = -nr/2, r? = (n - l)/2 - n T . (38) 
This gives a simple and very useful relation for Aq/A 2 ^, 

A 2 

-| = -0.5n T [1 - 0.5n T + 1.0(n - 1)] . (39) 
Substituting into Eq. (|36j ), we find the second-order consistency relation 

n T = -I [l-0.8n T + 0.15(n-l)], (40) 

or 

- = -7n T [1 + 0.8nr - 0.15(n - 1)] . (41) 

To the required order we can use the first-order truncation = — inside the brackets, 
thereby obtaining an alternative form, 



IT 

HT = ~7S 



T 

0.11— + 0.15(n 



(42) 



where nx is given in terms of the more accessible quantities (n — 1) and ^. 

Independent measurements of n, rix and ^ provide a powerful test of the inflationary 
hypothesis; in the space of these parameters inflationary models must lie on the surface 
defined by Eq. (14~l|) . In Figure 1 we illustrate the inflationary surface both without and 
with second-order corrections. The second-order corrections break the degeneracy in the 
(n — 1) direction, as well as typically reducing |J viewed as a function of n T and (n — 1). 
However, the portions of the surface that feature large corrections are not favored by present 
cosmological data, and further, are susceptible to higher-order corrections. (Indeed, well 
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away from scale-invariance the surface would be noticeably different even just using Eq. ( f42| 
instead of Eq. (f|l]), which differ by third and higher order terms.) 

Obtaining the reconstruction equations is simply a matter of substituting into Eqs. (3.4) 



(3.6) and (3.15) of Ref. []13[ for V, V and V" respectively. We give two alternative forms 
for each, the first using ut and the second substituting ^ for ut using the second-order 
consistency equation. They are 



1.65(1 - 1.4n T )T, 

1.65 (i + 0.2o|)t, 

±8.3^/^n^[l - l.ln T 
T 



0.03(n - 1)]T, 




1 + 0.21- -0.04(n- 1) 
o 



(43) 



(44) 



V&fm 



21 



in 



1) - 3n n 



lAn 2 T 



+0.6n T (n - 1) - 0.2(n - l) 2 + 1.1 



dn 
dink 



T. 



21 



n 



T 

0.43 ? 



° 073 'I)' 



-0.015^(n 



1) - 0.2(n - lY + 1.1 



dn 
dink 



T. 



(45) 



These expressions are accurate to second-order. Naturally, they agree with the first-order 
expressions given earlier. 

Though no expression is given in Ref. |13|] for V", by using the lowest-order expressions 
for e, n, and a third slow-roll parameter £, and Eq. (3. 13) which relates the three to dn/d In k, 
one can obtain the first-order expression, 



V"/m Pl = ±104v^ 



dn/dln k 



±104a 



'IT 
7S 



-7 



nx 

dn/dlnk 



6n T + 4(n - 1) 



T 



+ 0.9- + 4(,- 1 ) 



(46) 



where the overall sign is to be the same as that of V. The second-order term would require 
yet another observable. As remarked in Ref. |13|, even this first-order expression features 



the rate of change of the scalar spectral index, which is likely to be very difficult to measure. 
Realistically then, in the near term only the value of the potential and its first two derivatives 
are likely to be accessible to accurate determination. 

The final step in reconstructing the potential is to use d<j)/dN to the desired order, to 
find the range of that corresponds to the eight or so e-foldings of inflation relevant for 
astrophysics. To proceed, we may simply carry out a Taylor expansion of <p about 5O , to 
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whatever order we believe is appropriate, 



- 5 o =(N- 50) 



dN 



+ ■ 



(47) 



050 



This is a double expansion, in the sense that the coefficients are themselves obtained as a 
series expansion in the slow-roll parameters. 

To proceed, we use as a starting point the exact formula 



m 



PI TTl 



(4£ 



which, along with dN/dt 
be calculated, 



-H, yields the relation from which the Taylor coefficients may 



m P1 H' 



dN 4tt H' (49) 
To get a given coefficient in the Taylor expansion for (p^, one simply calculates d l <p/dN l 
expanding to the desired order in the deviation from scale invariance.0 For example, taking 
only the first term in the 0at expansion and working to first-order yields the expression 
already given in Section 1.1. We give the first coefficient in the <pN expansion to second- 
order and the second coefficient in the <pN expansion to first-order only, 



mpi 



>N - <?50 



± 



mpi 



V^nr [1 + 0.1n T + 0.1(n - 1)] (N - 50) 



4V87T 



yf^r [(n - 1) - n T ) (N - 50) 2 + 



(50) 



with both signs again agreeing with that of V. 

In the process of reconstruction, we shall use the first-order expansion for <Pn — 4>5o in 
first-order reconstruction, and the second-order expansion in second-order reconstruction. 



4 Reconstructing an exponential potential 

A useful testing ground for reconstruction is the exponential potential, the only known case 
where the perturbation spectra can be derived exactly analytically || [Tj|. For the potential 

VW=V a ev(-M-±), (51) 

V V V m PlJ 

6 Note this procedure differs slightly from that in Ref. fl^] , where dip/dN was expanded linearly about 
(j>50 and <Pn was solved for exactly, cf. Eq. (8). This results in an exponential, whose expansion picks up 
the (N — 50) and (N — 50) 2 terms correctly to lowest order in the deviation from scale invariance, though 
not the higher-order terms in the (N — 50) term which would require higher-order terms in the expansion of 
d(j)/dN . There is an overall sign error in Eq. (8) of Ref. fll^] . 
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the scale factor grows exactly as t p . Compared with the lowest order expressions, the ampli- 
tudes A and A T , or A 2 S and Aq, are both multiplied by the same p-dependent factor R 2 (p), 
where 

R(P) = 2^-» r|3/2 + 3 1 / / 2 ( ] p - 1 " (1 - , (52) 

where T(- • •) is the usual gamma function. Both scalar and tensor spectra are exact power 
laws with spectral indices (n— 1) = n T = —2/ (p— 1). The scalar- field solution is characterized 
by 

^ mP1 - # mP1 - V(^) = ^exp[2(JV-50)/p]. (53) 



4tt t diV v/llrp 
The expressions for T and S* can be obtained exactly by integrating Eqs. (^) and ([]), 

5 = 2.2/(, l )E 2 (p)> T , (54) 
m Pl x 50 

T/true 

T = O.ei^K)^ 2 ^)^, (55) 

where the numerical factors f(n) = 1 + 1.15(n — 1) + ■ • • and (?(^t) = 1 + 1.3nr + • • ■ arise 
from the n, tit dependence of the Sachs- Wolfe integrals, cf. Eqs. (||, H). 

We are now ready to carry out an array of reconstruction methods. Because we are 
using exact expressions to generate the spectra, this procedure is more ambitious, and more 
realistic, than those attempted thus far || [12| , where the trial spectra were produced using 



the slow-roll approximation. For the general inflationary potential, exact results are not 
known, and so this procedure is not possible^. However, our method here should give a more 
realistic estimate of inherent errors even in the general case. 

There are two distinct types of error. The first is error in the value of the potential at 
050, due to third-order and higher terms. By substituting the expression for T in Eq. flSo] ) 
into Eq. (|H]) or (|3|) for V50 we can compute that error: 

= 0(nr)(l - lAn T )R(p) 2 . (56) 

The second error involves the shape of the potential, which depends on the ability of the 
chosen expansion to match the potential over the eight interesting e-foldings. 

We have chosen as a specific example an exponential potential with p = 43/3. We did 
so because this leads to about the largest departure from scale invariance that can still be 
regarded as observationally viable, (n — 1) = nr = —0.15 and % ~ 1, and thus realistically 
represents the most challenging example of reconstruction. The exact potential is shown in 
Fig. 2 along with the results of five different reconstructions. 

To begin, consider the error in estimating V^ ne ; we have g(riT = —0.15) = 0.824 and so 
Vso/V^™ ~ 0.95, a modest 5% error due to the neglected higher-order terms. As we always 

7 Of course, we are going to pretend that we don't know the potential is exponential to demonstrate our 
methods. 
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include the second-order term in V50, the error is the same in every method we look at. Had 
the first-order expression for V 50 been used instead, corresponding to the neglect of the factor 
of (1 — IAut) in Eq. (|TDD, then the underestimation would have been about 20%. 

Let us now consider the shape, which we note depends on T and S only through their 
ratio. The important distinction between different methods is the difference in required 
input data; methods needing only n and ^ have the advantage of depending only on the 
information that is easiest to obtain. Requiring dn/dlnk in addition, while offering more 
accuracy, is setting a much trickier observational task, though upper bounds are also useful 
in the absence of actual determinations. 

As a starting point, let us take the equations derived in Ref. ||12|| , which are primarily 
first-order though they include the second-order correction to V50, cf. Eqs. ([T0HT2|). In this 
extreme example, the quadratic Taylor series based upon this does a bad job of approximating 
the shape of the potential, as it turns upward for large (0 — 5O ) due to the truncation at 
the (0 — 050 ) 2 term (see Fig. 2). 

If we now require knowledge of dn/dlnk, the Taylor series approach can be improved in 
two ways. We can now take V50, V§ , and to second-order; however, the improvement 
is rather minimal. Alternatively, we can stick to first-order expressions, but include the V§q 
cubic term. Again the improvement is modest, though at least the unwanted minimum has 
been eliminated. One could go further and take V50, V§ Q , and V$ Q to second-order and to 
lowest order, which we haven't illustrated, again seeing only modest gains for the increased 
observational requirement. 

The Taylor series having been unimpressive, let us progress in a different direction. With 
only n and ^ , as an alternative to the Taylor series one can construct the Pade approximant 
based upon it, taking V50 to second-order and V^ and V 5 ' to first-order. This represents 
a substantial gain on the Taylor series to that order without requiring any additional input 
information. With this minimal information, it is a much better method. Reintroducing 
dn/dlnk allows this method to be extended to second-order, where the reproduction of the 
shape of the potential is excellent. To include the third derivative term would necessitate a 
more complicated (non-diagonal) Pade approximant, which doesn't seem warranted at the 
moment. 

What is the upshot of this comparison? Recalling that we have chosen an example with 
extreme deviation from scale-invariance, the second-order corrections are reassuringly small 
and only improve the shape of the reconstructed potential slightly. The addition of the 
third derivative term in the Taylor series gives a slightly more significant improvement, but 
at the price of its dependence upon dn/dlnk even at lowest order. The most remarkable 
improvement involves the use of Pade approximants. Even without knowledge of dn/dlnk 
the shape of the potential is reproduced far better than with the higher-order Taylor series 
which does require that knowledge. As noted previously, the improvement results from the 
fact that the Pade approximant is not truncated; further, even in situations where truncation 
of the Taylor series does not lead to problems, the Pade approximant still proves valuable 
as its Taylor expansion coincides with that of the original expansion. We therefore conclude 
that Pade approximants provide a significant improvement in the perturbative reconstruction 
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of the inflationary potential. 



5 Discussion 

By presenting the second-order reconstruction equations directly in terms of observables, we 
have been able to assemble and to compare an array of different perturbative reconstruction 
techniques based upon cosmological observables. Our work extends previous work in several 
important ways. 

First, we have placed the perturbative reconstruction process on a firmer foundation by 
addressing the important issues of convergence and term ordering. We have emphasized 
that the observational data themselves can be used to decide whether or not perturbative 
reconstruction is well justified and sensible. In particular, we have shown that the Taylor 
series for the potential is absolutely convergent and that terms in the expansions for the 
observables must decrease in size as the number of derivatives increase for the case where 
the spectral indices do not vary significantly over the astrophysically interesting scales, or, 
if they do, their absolute change is small. 

Perhaps our most interesting result is the introduction of the Pade approximant as an 
alternative to the Taylor series in perturbative reconstruction. It can be obtained from a 
Taylor series regardless of the order (in the deviation from scale invariance) to which the co- 
efficients of the Taylor series has been obtained. In our worked example, the improvement in 
reproducing the shape of the potential as compared to the Taylor series is striking, especially 
considering that no extra observables are required. 

We have shown that the second-order corrections to the Taylor series coefficients are 
generally small, and that those for V50 and V 5 ' only depend upon the same quantities as 
the first-order expressions (S, T, and n). The corrections to V 5 'o however require a new 
observable such as dn/dlnk, and by deriving for the first time an explicit expression we have 
confirmed that even the lowest-order term in V 5 'q requires this challenging observable. 

Finally, one of the most important aspects of reconstruction is that it is overdetermined: 
Any set of cosmological observables supplies degenerate information regarding the potential 
and its derivatives, thereby providing an important consistency check. In particular, the 
tensor spectral index can be expressed to second-order in terms of S, T, and n by the 
relation: tit = — \ §[1 + 0.11^ + 0.15(n — 1)]. In cases that are observationally viable, the 
second-order corrections are small. 
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Appendix: Some Relations between Notation 

For the convenience of the reader, we summarize here some relations between the notation 



used here and that in Ref. fL3 |, from which several important results were taken. In that 
paper, the spectra As and Aq were defined so as to include any scale-dependence within 
them, i.e., they are functions of k. In circumstances where the spectra can be approximated 
by power-laws, these are related to the amplitudes A and At in this paper, which are just 
numbers, by 

Aik/ho)- 1 = ^A&k), (57) 



A T (k/hoT T = 2A%{k). (58) 

Even in cases where the spectra cannot be described by power-laws, the correspondence 
holds at k — kso. 



In Ref. |L3| , slow- roll parameters e and rj are introduced, 

which are again in general fc-dependent. As indicated in Section 2 of the present paper, they 
can be related to the spectral indices to various orders, e and r\ being of the same order in 
perturbation theory as (n — 1) and tit- To lowest-order they are constant, corresponding 
to power-law spectra. At lowest-order e = 16ttx 2 , but higher order corrections break this 
relation. 
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Figure Captions 

Figure 1: The consistency plane for inflation in n — ut — ? space, the flat surface being the 
lowest-order result and the curved one incorporating the second-order corrections, given by 
Eq. (HD. 

Figure 2: An array of different reconstructions of an exponential potential with (n — 1) = 
ht = —0.15 (p = 43/3). The longer dotted line indicates the exact potential. The three 
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different line styles correspond to three different reconstruction strategies; solid is Taylor 
series truncated at (0 — 05o) 2 > dashed is Taylor series truncated at (0 — 05o) 3 and dash- 
dotted is the Pade approximant based on the former of these. The upper line of a given 
style uses coefficients to first-order in the deviation from scale invariance (save V50, which is 
always second-order), while the lower, where plotted, is second-order in all coefficients. The 
length of the curves corresponds to eight e-foldings. 
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